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^ The SU(2) gauge theory with two massless Dirac flavors constitutes the building block of several 

models of Technicolor. Furthermore it has also been used as a template for the construction of a 

^ natural light asymmetric, or mixed type, dark matter candidate. We use explicit lattice simulations 
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to confirm the pattern of chiral symmetry breaking by determining the Goldstone spectrum and 
therefore show that the dark matter candidate can, de facto, be constituted by a complex Goldstone 
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^ spectrum and demonstrate that it is well separated from the Goldstone spectrum. 
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I. INTRODUCTION 



The nature of dark matter (DM) is an important problem in modern physics. DM plays 
a key role in the formation of large structures and the evolution of the Universe. It is 
also widely expected to provide a link to physics beyond the Standard Model (SM). For 
these reasons much experimental, observational, and theoretical effort has been devoted 
to shedding light on DM. It is popular to identify DM with Weakly Interacting Massive 
Particles (WIMPs). 

Many properties of a WIMP are not constrained by our current knowledge of DM, for 
example the WIMP may or may not be a stable particle [lj and it may or may not be 
identified with its antiparticle [2]. Asymmetric DM refers to a scenario where the WIMP's 
antiparticle has been annihilated away, leaving the WIMP itself as the observed DM. 
Asymmetric DM candidates were put forward in [2] as technibaryons, in Q as Goldstone 
bosons, and subsequently in many diverse forms ||4l-[10|. There is also the possibility 
of mixed DM |[TT|, i.e. having both a thermally-produced symmetric component and an 
asymmetric one. 

Null results from several experiments, such as CDMS |fl2fl and XenonlO/100 |fT3l 



[14| , have placed stringent constraints on WIMP-nucleon cross sections. Interestingly 
DAMA fT5| and CoGeNT ||T6l have both produced evidence for an annual modulation 
signature for DM, as expected due to the relative motion of the Earth with respect to the 
DM halo. These results support a light WIMP with mass of order a few GeV, which offers 
the attractive possibility of a common mechanism for baryogenesis and DM production. 
At first glance it seems that the WIMP-nucleon cross sections required by DAMA and 
CoGeNT have been excluded by CDMS and Xenon upon assuming spin-independent 
interactions between WIMPs and nuclei (with protons and neutrons coupling similarly 
to WIMPs), but a number of resolutions for this puzzle have been proposed in the litera- 
ture |fl7l - l22| . Interestingly, new results from the CRESST-II experiment report signals of 
light DM ||Z3]|. 

A composite origin of DM is an intriguing possibility given that the bright side of 
the universe, constituted mostly by nucleons, is also composite. Thus a new strongly- 
coupled theory could be at the heart of DM. In this work we investigate, for the first 
time on the lattice, a technicolor-type extension of the SM expected to naturally yield a 



light DM candidate, as introduced in [8] and used in [|22| to reconcile the experimental 
observations. 

II. A LIGHT DARK MATTER TEMPLATE: SU(2) TECHNICOLOR WITH TWO DIRAC FLA- 
VORS 

Asymmetric DM requires at least one complex field with a (nearly) conserved Abelian 
quantum number. It would be exciting if the theory predicting the existence of this state 
is also directly involved with the breaking of electroweak symmetry. Therefore a natural 
candidate is a technicolor model. Most of the states of the theory are much heavier than 
a few GeV, but if there are Goldstone bosons not absorbed by the longitudinal degrees of 
freedom of the SM massive gauge bosons, they become primary candidates for producing 
a natural hierarchy between the GeV and the TeV scale. 

The minimal technicolor theory that breaks the electroweak symmetry and features 
a light DM state was constructed first in [|24ll25B . This theory is SU(2) technicolor with 
two Dirac flavors, which has global symmetry SLT(4) expected to break to Sp(4). Five 
Goldstone bosons are generated of which three are eaten by the SM gauge bosons and a 
complex one (a techni-diquark which is essentially a technibaryon of this SU(2) theory) 
remains in the spectrum. This is our candidate for the light asymmetric DM particle. 

The walking version of this model, known as Ultra Minimal Walking Technicolor 
(UMT), has been constructed in [8] and it contains, besides the fermions in the funda- 
mental representation, also a Dirac fermion in the adjoint. For the present exploratory 
work, we begin with Section III of flSJ but without the adjoint fermion. Moreover, when 
implementing our lattice simulations we make two additional modifications: all elec- 
troweak interactions are omitted from the lattice simulations because that physics is 
well-understood by perturbative methods, and explicit technifermion masses are added 
because lattice simulations with exactly massless fermions are not practical. The resulting 
Lagrangian is very simple 

X = + U(iy i D, - m)U + D(i 7 ^D, - m)D , (1) 

where U and D are the two techniquark fields having a common bare mass m, F" is the 
field strength, and is the covariant derivative. The Dirac and technicolor indices of U 

3 



and D are not shown explicitly. 

Lattice simulations will be used to extrapolate to m = 0, and in that limit the Lagrangian 
has a global S17(4) symmetry corresponding to the four chiral fermion fields 

U L = ^(l-y 5 )U, U R = ^(l + y 5 )U, D L = l -{l-f)D, D R = ^(l + y 5 )D. (2) 

For m + 0, the St7(4) symmetry is explicitly broken to a remaining Sp(4) subgroup as 
follows. The Lagrangian from ([T| can be rewritten as 

X = ~F; IV F^' + iUy^U + iDytDpD + ^Q T (-io 2 )C EQ+^ (Q T (-io 2 )C EQ) + , (3) 

where 
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C is the charge conjugation operator acting on Dirac indices, and the Pauli structure -io 2 
is the standard antisymmetric tensor acting on color indices. Under an infinitesimal SLT(4) 
transformation defined by 



the Lagrangian ^ becomes 



Q 
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Q. 



(5) 



X -> X + y Yj a "Q T (- i(j2 ) c { ET " + T " TE ) Q + h - c - / 



(6) 



n=l 



where T n denotes the 15 generators of SLT(4) and a n is a set of 15 constants. The only 
generators that leave the Lagrangian invariant are those that obey 



ET n + T nT E = 



(7) 



which is precisely the definition of an Sp(4) Lie algebra. From this, it is straightforward 
to derive the ten Sp(4) generators in a specific representation; see the appendix of [8J. 

For m = the Lagrangian retains the full SU(4) symmetry but, by analogy with the 
SU(3) theory of QCD, one might expect dynamical symmetry breaking associated with 
the appearance of a nonzero vacuum expectation value, 



(1717 + DD) + 
4 



(8) 



Since this vacuum expectation value has the same structure as the terms containing m 
in the Lagrangian, the dynamical breaking would also be SLT(4) — > Sp(4). According to 
Noether's theorem, the five broken generators would be accompanied by five Goldstone 
bosons. 

Of course this suggestion of dynamical symmetry breaking must be checked nonper- 
turbatively using first-principles lattice simulations. As reported below, we have done 
so and our lattice simulations provide direct verification of this dynamical symmetry 
breaking. 

There have been several previous lattice studies of SU(2) gauge theory with fermions 
in the fundamental representation [|26|, mainly motivated by interest at nonzero chemical 
potential, but all of these studies relied on the staggered action where the number of 
fermions is a multiple of 4. Our minimal technicolor theory requires Nf = 2 and thus 
we use the Wilson action instead of staggered fermions. When studying chiral symmetry 
breaking scenarios with Wilson fermions attention must be paid to the presence, on a 
lattice, of the unphysical Aoki phase. For fixed gauge coupling, the Aoki phase is entered 
as the quark mass is reduced. An analytic discussion of the Aoki phase symmetries for 
this theory is provided in [273, and three lattice studies are also available ||28] - l30| . For our 
present simulations, we avoid the Aoki phase and work exclusively in the physical phase. 
There have been very few previous lattice results reported for SU(2) gauge theory with 
2 fundamental fermions II30H32M and in each case the primary focus was on a different 
action (either Nf > 2 or adjoint representation fermions). Our work is the first lattice study 
focused on the mass spectrum of the two-color two-flavor theory, which is the familiar 
technicolor template. 

III. LATTICE HADRON OPERATORS AND EFFECTIVE FIELD THEORY 

The creation operator for a meson is the Hermitian conjugate of its annihilation oper- 
ator, and a set of local annihilation operators for mesons is 



Og = U(x)TD(x) , 
Og, = D(x)TU(x) , 




(9) 
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where T is a chosen Dirac structure. In this work we consider T = 1, y 5 , y^, or y^y 5 . In 
lattice simulations, meson masses are extracted from the time dependence of correlation 
functions, for example 



x y 

= E E Tr [my)D(x)y°TYU(x)U(y)] , (10) 

* y 

where Tr[- • • ] denotes a trace over Dirac and color indices and we have dropped the 



vacuum expectation values for the propagators as explained in Appendix A 

Perhaps surprisingly, local diquark (i.e. baryon) correlation functions provide no new 
data in this theory. To understand why, notice that the available local diquark operators 
are 

OB, = U T (x)(-io 2 )CTD(x) , 



UD 

Og, = D T (x)(-ia 2 )CTU(x) . (11) 



(Some operators containing U T ■ ■ ■ U or D T • • • D are identically zero.) The diquark corre- 
lation function is therefore 



- h) = E E °uoM K>)) + 

x y 

(12) 



£ £ Tr [rD(y)D(x)/r + c + (-w 2 )V 0T w T WLr T (y)(-w 2 )c 



x y 



and this can be rewritten by using two properties of the charge conjugation operator: one 
for a Dirac matrix, 

y lT = -Cy f 'C + , (13) 



and the other for the Wilson fermion matrix, 

. T 



to arrive at 



C- 1 (-w 2 )- 1 (LT(y)LZ(x)j C(-io 2 ) = U(x)U(y) , (14) 



C ud^ ~ h) = E E Tr [my)D(x)y°TyU(x)U(yy] 

x y 

= Cg(tx ~ ty) (15) 



for any choice of T. This conclusion means that a lattice simulation will find numerically- 
identical correlation functions, and thus identical masses, for the meson and diquark. We 
have verified this explicitly in the lattice simulations described below. 

Notice that the degenerate pairs have equal angular momentum but opposite parities, 
for example 

/MS,) =/(<>£). a*) 

P(Og,) = -P (<*>). (17) 

This relationship for } p has also been mentioned, for example, in ||33| . We must conclude 
that the three pseudoscalar meson Goldstone bosons are accompanied by two scalar diquark 
Goldstone bosons. This is in contrast to the identification made in [8], where it was 
assumed that all five Goldstones would be pseudoscalars. Despite this, the effective 
Lagrangian in [8] remains unaltered because all scalar and pseudoscalar particles were 
retained in that work with the correct assignment with respect to the broken and unbroken 
generators of the chiral symmetry groupQ 

Now that the five Goldstones have been identified, we can define an effective theory 
wherein those five are the only fields. Indeed, our lattice simulations (discussed below) 
indicate that the non-Goldstone scalar/pseudoscalar hadrons are even heavier than the 
vector mesons, making it quite natural to integrate all non-Goldstones out of the effective 
theory. The annihilation operators are conveniently collected into a compact form, 

6£ = L (Q T Htf 2 ) crr q)° (i> . (is) 

where O^" is a generic name for the n'th particle with Dirac structure T. This SX represents 
the effective Lagrangian couplings for external fields 0^ r) ". For the special case of the five 
Goldstone bosons, we sum n only over the five broken generators and we must choose 
T = y 5 . The result is 

5£g = Q\-io 2 )Cy 5 QQ (19) 
where the matrix Q contains the five Goldstone bosons (i.e. the pseudoscalar mesons n + , 

1 In practice the only change is the redefinition of YhjD and a * the underlying operator level with YhjD 
and TItjp in |8J. 
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IT , n° and the scalar di quarks Hud, TLjp): 






V2n UD 


n° 


V2rr 


-V2n UD 





V2rr 


-n° 


-n° 


- V2rr 







- V2n + 


n° 


V2n^ 






(20) 



The matrix Q is closed under Sp(4) transformations: the Goldstone bosons form a five- 
dimensional representation of Sp(4). The matrix Q is also identical to M4 as defined by 
equation (22) of [8 J after the non-Goldstone fields are removed from M4. 

To summarize, an effective field theory for the five Goldstone bosons is obtained by 
using the matrix Q in place of M4 in J8J. The three pseudoscalar Goldstones (IT*, n°) 
are responsible for electroweak symmetry breaking, and the two scalar Goldstones (Hud, 
TL-jjp) are the DM candidate and its antiparticle. 



IV. METHODS FOR NUMERICAL SIMULATIONS 

The standard Wilson action, 
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Sw = 2 Tu f 1 ~ 5 ReTrL ^ (*) u v(* + + v)Ul (x)) + (4 + mo) 

y^U^x^Cx + A) + ${x + + y f ,)*i>W*)l (21) 



is used for this study of SU(2) gauge theory with two mass-degenerate fermions in the 
fundamental representation. Configurations were generated using the HMC algorithm as 



implemented in the HiRep code ||34f. A total of 12 ensembles were created, corresponding 
to six different bare quark masses m for each of two different gauge couplings jS; see Table 
[I] for details. We note that B30l also contains some simulations of this theory at f> = 2.0 
and our findings for the pseudoscalar meson mass and PCAC quark mass are consistent 
with that paper. For other parameter choices, see B3TI and II32II . In the present study, 
all lattices are L 3 xT = 16 3 x 32 with periodic boundary conditions in each direction. 
Every ensemble contains 35 configurations separated by 20 unused configurations after 
an initial thermalization of 320 configurations. Figures [T] and [2] display the fluctuations in 
the average plaquette within selected Markov chains. 



TABLE I: Numerical values of the gauge coupling and the bare quark mass used to 
generate the 12 ensembles of this project. 





m 








2.0 


-0.85, 


-0.90, 


-0.94, 


-0.945, -0.947, -0.949 


2.2 


-0.60, 


-0.65, 


-0.68, 


-0.70, -0.72, -0.75 



0.575 




configuration number from an ordered start 

FIG. 1: Evolution of the average plaquette for selected bare masses at jS = 2.0. Our 
ensemble of 35 configurations is comprised of those numbered 320, 340, 360, . . . 1000. 

Quark propagators are created from wall sources built of random U(l) phases at ev- 
ery lattice site on one chosen time step, tj. Inversions are performed with the standard 
BiCGstab algorithm. To reduce statistical fluctuations, correlation functions are averaged 
over all lattice times < t t < 31. Correlation functions then depend on the time tt cor- 
responding to the annihilation operator. This annihilation operator is local and summed 
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0.62 




configuration number from an ordered start 

FIG. 2: Evolution of the average plaquette for selected bare masses at f> = 2.2. Our 
ensemble of 35 configurations is comprised of those numbered 320, 340, 360, . . . 1000. 

over all spatial sites to produce a zero-momentum hadron. 
Our multi-state fits, 

it 

C(f) = £ fly COSh 2))' ( 22 ) 

7=1 

use every nonzero time separation t = tt—ti > 0, which avoids the subjectivity of choosing 
a fitting window. Pseudoscalar, vector and temporal-axial (A4) correlators are fitted to 
three states, n = 3 in ( |22| >; scalar and spatial-axial (Ai,A 2 ,A 3 ) correlators use n = 2. 
Statistical uncertainties are produced from each 35-configuration ensemble by creating 
150 bootstrap ensembles (having 35 configurations each). 

Isovector hadrons are sufficient for most topics addressed in this work, but isoscalars 
are used to study a specific issue. Isoscalar operators require the computation of single- 
site propagators, i.e. quark propagators that begin and end on a single lattice site. In 
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FIG. 3: The t dependence from the right-hand side of (23 1 for the case of [5 = 2.2 and 
mo = -0.75. The horizontal line is a fit to Euclidean times between 10 and 22 inclusive. 



practice, to obtain a signal requires an average over many sites, so isoscalar studies are 
expensive. Improved methods have recently been developed ||35ll36| , but the needs of our 
present project are fulfilled by a simpler tactic. We choose our "most physical" ensemble 
(jS = 2.2 and ra = -0.75) and calculate 2048 single-site propagators, that is 64 single-site 
propagators spread uniformly across each lattice time step. 

Two additional quantities that can be derived from two-point correlation functions are 
also valuable for the present work. One is the PC AC quark mass defined by 

r ( (A 4 (t + 1)P(0)) - (A 4 (t - 1)P(0)) \ 

m < = & ( 40Wb» ) ' (23) 

and in practice we can average over time separations satisfying 12 < t < 20. Figure|3]shows 
the example of /3 = 2.2 and m = -0.75. The other additional quantity is the Goldstone 
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boson decay constant that we extract from a simultaneous fit to three correlators 113711 , 




where L = 16 and T = 32. As stated above, we use n = 3 and fit all time separations 
except t = 0. The fitting parameters are ntj, mj/(m q Z P ), and fj/Z A . Notice that such fits 
only provide the decay constant divided by its renormalization constant Z A . Since Z A will 
approach unity in the continuum limit, we will study the ratio fj/Z A rather than fj itself 
in this exploratory study 



V. RESULTS FROM NUMERICAL SIMULATIONS 



As familiar from QCD, the axial Ward-Takahashi identity incorporates a partially- 
conserved axial vector current (PCAC) which defines a renormalized quark mass in 
terms of the axial current's derivative. Use of the explicit definition in p3| produces the 
relationship between renormalized and bare quark masses as plotted in Fig. |4j Over our 
range of bare masses, the PCAC mass is seen to be effectively linear in the bare mass. 

Goldstone bosons should be massless when the renormalized quark mass is zero, and 
should approach that limit according to 

mn oc ^m~ q . (27) 

The lattice data confirm this behavior, as shown in Fig. |5j At /3 = 2.0 our entire range of 
masses satisfies d27k. At /3 = 2.2 our range of mass values is broad enough that the two 



heaviest quarks show curvature coming from corrections to (27 1, and our lightest quark 
displays a finite-volume correction. We wish to emphasize that these plots in Fig. [5] are 
obtained using any of the five operators from Q and (11). All five Goldstone bosons 
are exactly degenerate. No single-site propagators were needed for computing these 
correlations functions. 
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FIG. 4: The PCAC quark mass as a function of the bare quark mass for (a) j6 = 2.0 and 
(b) f> = 2.2. In each case, the line is a fit to all six data points. 



The physical scale of technicolor can be specified by requiring that the Goldstone decay 
constant matches the observed electroweak energy scale, 



fn * 246 GeV. 



(28) 



Comparison to a lattice determination of fn would thus provide a direct interpretation 
of the physical scale for each lattice simulation. Up to the renormalization constant Z A , 
we can produce these lattice results by using (24|26]), as displayed in Fig. [6} A curvature 
is clearly visible in these plots of fn/Z A versus m q , and for our purposes it is sufficient to 
notice that fnl^A ~ O(0. 1) in our simulations, i.e. an order of magnitude below the lattice 
cut-off. 

Figure |7j besides displaying again the Goldstone isovector spectrum, can be used to 
verify that non-Goldstone type states, such as the scalar and spin one isovector particles, 
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FIG. 5: The Goldstone boson mass squared as a function of the PCAC quark mass for 
(a) jS = 2.0 and (b) /3 = 2.2. In (a) the line is a one-parameter fit to all six data points, but 
in (b) the line is a fit to the three points within 0.05 <m q < 0.15. 



remain massive in the m q = limit. Much like QCD, the vector meson mass is nearly 
an order of magnitude larger than fij, but significantly lighter than the isovector scalar 
and axial vector mesons. Lattice artifacts are presumably sizable for masses m > 1, but 
our results do extend into the region where all masses are less than 1. Notice the exact 



degeneracy between mesons and diquarks of opposite parity as anticipated in Section III 



Although isoscalar correlation functions bring the major challenge of single-site prop- 
agators, we want to consider the specific case of the isoscalar pseudoscalar meson. Our 
goal is to verify explicitly that this meson is not an additional Goldstone particle. Figure [8] 
displays the results. Panel (a) is the correlation function for the isovector, i.e. the quark- 
connected diagram only A multi-exponential fit to those data produce the corresponding 
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FIG. 6: The Goldstone boson decay constant as a function of the PCAC quark mass for 
(a) j6 = 2.0 and (b) j8 = 2.2. 



data point in Fig. |5j Panel (b) is the contribution from the quark-disconnected diagram 
built from a pair of single-site quark propagators. Panel (c) is obtained from the upper 
two panels according to // (c)=(a)-2(b)", where the factor of 2 accounts for UU plus DD 
contributions and the minus sign accounts for anticommuting fermions. In panel (c), data 
at all timesteps between 7 and 24 are statistically consistent with zero, and time steps 
near the source at t = correspond to states much heavier than the Goldstones. Thus 
we see how the quark-disconnected diagrams have canceled the Goldstone signal out of 
the isoscalar. This is the expected conclusion since the symmetry breaking StT(4) — > Sp(4) 
requires exactly five Goldstone bosons. 
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FIG. 7: Mass spectrum of isovector particles, with various } p quantum numbers, as a 
function of the PCAC quark mass for (a) jS = 2.0 and (b) /3 = 2.2. 



VI. CONCLUSIONS 

From many previous studies pHTTj the lightest neutral technibaryon has become a 
classic prototype for asymmetric DM, and |[3j, HUB! E21 focused on the possibility of that 
technibaryon being a pseudo-Goldstone boson. In this work, we have performed lat- 
tice simulations of SU(2) gauge theory with two techniquarks which has been used as a 
template for the construction of composite light DM models ||8l 122]. For the first time, 
we verify the existence of the pseudo-Goldstone state to identify with the composite DM 
candidate. In fact, besides the standard three Goldstone bosons, to be eaten up by the 
SM gauge bosons, we observe two more Goldstone technibaryons. These latter states 
are scalars rather than pseudoscalars. The masses of other phenomenologically relevant 
technihadrons were found to be an order of magnitude heavier than the Goldstone decay 
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FIG. 8: The (a) connected, (b) disconnected, and (c) total contributions to the 
correlation function of an isoscalar pseudoscalar meson at /3 = 2.2 and mo = -0.75. 



constant, thus justifying the construction of an effective field theory containing only the 
five Goldstone fields introduced in [|24| and then adapted to UMT in [8J. We have estab- 
lished the expected pattern of chiral symmetry breaking in the absence of the electroweak 
interactions. Typically these interactions tend to destabilize the vacuum (see |38| for a 
recent study) which can be, however, re-stabilized when extending the technicolor theory 
to give masses to the SM fermions and to the otherwise massless DM candidate. 

One might ask whether a light isoscalar scalar meson is present in this gauge theory, 
reminiscent of the cr//o(600) meson of QCD, that could be viewed as a composite Higgs 
boson. The isoscalar scalar channel is a major challenge for lattice simulations B39l , and 
we must leave this question to be the focus of a future study. 

In the future it will be very interesting to determine more accurately the mass spectrum 
of the theory, but there are many additional topics waiting to be explored as well. Here we 
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mention three of them. The effects oifour-fermion interactions, relevant for a realistic exten- 
sion of the SM, should be investigated ||40ll . The extension to nonzero temperature would 
lead to a characterization of the early universe electroweak phase transition. Knowledge 
of its order and strength is essential for investigating the possibility of providing correct 
baryon and asymmetric DM genesis at the electroweak scale [|4T] - |43| . The nonzero matter 
density regime is also useful for asymmetric DM physics El |8j given that an asymmetry 
implies the presence of a nonzero chemical potential. In contrast to QCD, the extension 
to nonzero matter density is directly possible in lattice simulations of this SU(2) theory. 
That is because the fermions belong to a pseudoreal representation of the SU(2) gauge 
group and therefore a chemical potential does not lead to an imaginary action. At nonzero 
matter density a number of interesting phases can emerge, such as baryon superfluidity 
(due to the condensation of the diquark technibaryon) and rotational invariance breaking 
(due to spin-one condensates [44 -46] arising from Lorentz symmetry breaking at nonzero 
chemical potential). 
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Appendix A: Proof of Equation p4| 

The Wilson fermion matrix as obtained from Eq. ( |2T] > is 

(UiyMx))- 1 = (4 + m )5 xy - ~ £ ((1 - y,)U,(y)5 y+ ,, x + (1 + y^ZljMV^) . (Al) 

f 

The transpose is 

([(UiyMx)] 7 )- 1 = (4 + m )5 xy - i £ ((1 - yJ)ll£(y)S y+ ^ + (1 + yfiU^xfy^ . (A2) 

Applying the charge conjugation operator C we obtain 

C-HMyMx)] 7 )-^ = (4 + m )5 xy - ^ £ ((1 + y,)(iJ(y)S w + (1 - y,)U;(x)5 lJiXn )) . (A3) 

f 
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Recalling that for a general SU(2) matrix H the following identity holds: 

{-io 1 T 1 V{-io 2 ) = U , 

we deduce: 

C- 1 (-w 2 )- 1 <[LZ(y)H(x)] T )- 1 C(-zcj 2 ) 

= (4 + m )5 xy ~\Yj{ {1 + >V) U J(y)<W + (! - y^M'W) • (A4) 

f 

By direct comparison with Eq. ( |A1 > we arrive at 



C- i (-w 2 )- i <[LZ(y)U(x)] J )C(-^ 2 ) = (U(x)U(y)) , (A5) 



which corresponds to Eq. (14). In the main text we dropped the (..) symbols, indicating 



the vacuum expectation value, to ease the notation. 
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